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Effectiveness of using laser field to produce entanglement between two dipole-interacting identical 
two-level atoms is considered in detail. The entanglement is achieved by driving the system with a 
carefully designed laser pulse transferring the system's population to one of the maximally entangled 
Dicke states in a way analogous to population inversion by a resonant 7r-pulse in a two-level atom. It 
is shown that for the optimally chosen pulse frequency, power and duration, the fidelity of generating 
a maximally entangled state approaches unity as the distance between the atoms goes to zero. 

"dv 

0^ , With recent experimental advances in the methods for coherent manipulation of quantum system on the level of 
individual particles, many of previously purely speculative problems become surprisingly up-to-date. In particular, 
much activity of physicists from different research fields is currently devoted to clarification of the entanglement 
concept ways for its quantification, purification, and creation. Casual creation of entangled states of atoms by 
coherent manipulation with light currently poses one of the biggest challenges in the field of quantum optics pj. 
^ | , Conversely, the resonant dipole-dipole interaction (RDDI) and cooperative relaxation effects associated with it are 
rather traditional topics of research. Recently the RDDI has been investigated as a source of interference phenomena 
in emission spectra |j,^|,|5| , super- and sub-radiance , photon bunching [|| , collisions in the laser cooling processes 
C*~) ' anc ^ as a mecnan i sm f° r realizing quantum logical gates pj] , p^ |). In this paper we address the question: how 

J> ' effectively can the RDDI along with coherent laser pulses be used for creation of multi- atomic entangled states. 
\Q '• 

In our model two identical two-level atoms are located at a fixed distance R and can be driven by laser beam that 
is either parallel or perpendicular (according to geometries identified, respectfully, as antisymmetric and symmetric, 
Fig. 1) to the radius vector R connecting the atoms. Within the interaction picture and rotating wave approximation, 
the evolution is described by the following master equation [^J : 
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^ . is the effective Hamiltonian describing the atoms' self-evolution and interaction with the laser field. Here 5 = ujl ~ uj a 
is the laser detuning from the atomic transition frequency, Qi,% are the complex laser driving Rabi frequencies for 
each atom, of, of, of, a i — of ± iaf, i — 1,2, are (using the well-known analogy between two- level atoms and 



spins) the spin-i Cartesian component and transition operators of the z-th atom, and we define g, / and 7 so that 
7n = 722 = 7, 712 = 721 = 37, X = /T< The distance dependent parameters g and /, describing, respectfully, 
the photon exchange rate and coupling due to the RDDI, are defined differently for different types of the atomic 
transition in question. Defining p = 0, q = 2 for Am = and p = 1, q = —1 for Am = ±1 transitions (with the 
quantization axis coinciding with R) and assuming that dipole matrix elements of the atoms are collinear to each 
other and perpendicular to R, we have the following expressions for g(<p) and /(</?) Jl3| | : 
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where ip = kR is the "dimensionless distance" between the atoms and k — uo a /c. Throughout the rest of the article 
we will consider Am — ±1 case for determinacy as for the Am = case the results are qualitatively same. 

It can easily be shown that the Dicke states, \ip e ) — \e) l |e} 2 , \ip g ) — \g) 1 \g) 2 , \ip s ) — -±=(\g) 1 \e) 2 + |e} 1 |3) 2 ), and 

\ip a ) = -i=(|(j() 1 |e} 2 — \e) 1 \g) 2 ) (where | e) ^ and \g) i are the upper and lower levels of the zth atom) are the eigenvectors 

of TL c s (when excluding the laser driving) , while the rest of the atomic dynamics can be described as radiative decay 
to/from the antisymmetric \ip a ) state with the rate j— = (1 — 3)7 and to/from the symmetric \ip s ) one with the rate 
7+ = (1 + 9)1 ■ 
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Analytical stationary solutions of ([!]) can be found for the case of symmetric excitation fii = ^2 = ^ (without loss 
of generality we can assume that Q here is real and positive) with the stationary populations of the Dicke states given 
by 
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Graph N S (Q,S), corresponding to ip = 0.5, is shown in Fig. 2a. The antisymmetric case, when the laser beam is 
parallel to i?,, allows no such simple analytical solution since in this case the relation between the Rabi frequencies 
for two atoms is more complex f2i = e l T22 = f2 (although Q here is again real and positive). However, the numerical 
solution is easily obtained and the corresponding dependence N a (fl, 6) is shown in Fig. 2b. 

If we aim to transfer the maximum amount of population into one of the maximally entangled states \^B a ) or 
|\& s ) by a short coherent pulse, a good criterion for finding optimal values of the laser field parameters, Q and 5, is 
whether the population of the corresponding level is close to 0.5 in the stationary solution. From the analysis of (Q) 
and the graphs in Fig. 2 we deduce that the optimal parameters can be well approximated by S op t = ±x((^)/2 and 
fl op t — \/\x{ l P)\l± with the upper/lower sign for the symmetric/antisymmetric geometries, respectfully. We then 
have \S opt \ — |xl/2 3> £l pt ^ 7± f° r sufficiently small distances, so that the transition of interest is saturated while 
at the same time the Rabi frequency is moderate enough to avoid broadband excitation of the Dicke level we are not 
interested in. 

Given the optimal parameters obtained in the previous section we proceed to find the fidelity of creation of the 
maximally entangled states, that is the maximum amount of population one can transfer using pulses of radiation. 
Considering the dynamics of the populations under optimal parameters laser driving that is turned on at the time 
instant t — (when all of the population is in \^ g ) state), we define the optimal pulse duration as the time when the 
population of the state we are interested in, reaches its first maximum. For the so chosen parameters we plot in Fig. 
3a the populations achieved by applying the optimal pulse, as a function of interatomic distance tp. The populations 
approach unity for both geometries as ip goes to zero suggesting that almost perfect transfer of population is attainable 
at small interatomic distances. 

It is easy to explain this result. As R goes to zero the energy splitting between the symmetric and antisymmetric 
state, equal to 2\ — 2/7, grows to infinity (of course, within the evident limitations of negligibility of exchange effects 
M ) . The "parasitic" excitation of the levels we are not interested in is then avoided by shifting the laser frequency so 
that we are in resonance only with one of the excited Dicke states. And if we have in possession arbitrarily strong and 
arbitrarily tunable laser, as ip goes to zero we can produce shorter and shorter pulses thereby decreasing the decay 
probability during the pulse. 

In addition, in the case of the symmetric excitation laser driving matrix element for the transitions involving the 
antisymmetric state vanishes. This means that its population can only come from decay of the upper J\I/ e ) level. 
In fact, in the stationary solution (Q) the two states even have the same populations, which are negligible for large 
detunings. In the antisymmetric excitation case, however, the situation works against us. As the interatomic distance 
goes to zero the Rabi frequencies of the two atoms become closer in phase, diminishing along the way the matrix 
elements of transitions involving antisymmetric state to zero. But even with such a small value of excitation efficiency, 
we can still manage to transfer the population to the antisymmetric state because the symmetric excitation remains 
far off resonance, and get as a reward the increased lifetime r — 1/ ((1 — g(R))~/) 3> 1/7 of our maximally entangled 
state (using this "durability" of the antisymmetric \^ a ) state, we can also create some entanglement passively as 
described in E3, but it is difficult to obtain a high fidelity that way). 

Now, instead of calculating different entanglement measures |l5|] (and then figuring out which of them better suits 
our purposes), we will consider how much our created states violate a simple Bell inequality. It can be easily shown 
jjj that for any classical local variable distribution the probabilities of finding the atoms' "spins" aligned along n z 
direction after coherent rotations fulfill the following Bell inequality: 

P diff (0, 2tt/3) + Pdiff (2vr/3, -2tt/3) + P diff (0, -2tt/3) > 1, (5) 

where the Pdis(<pi, ^2) is the probability of getting different results of the measurements of the two spins, i.e., of 
finding one spin aligned along the n z direction and the second one against it after the first spin is rotated by <p\ and 
the second one by ip2 around OX or OY axis. In quantum mechanics, however, the left hand side of d|) amounts 



2 



to only 0.75 for the pure maximally entangled \^ a ) state. To apply the same Bell inequality (|J) to the case of the 
symmetric excitation geometry we first perform a 180° rotation along n z -axis with one of the atoms (thus turning 
\^ s ) state into \^ a ) one) and then perform the measurement as described above. In Fig. 3b the l.h.s. of (||) is plotted 
against the dimensionless distance ip showing that for <p less than rs 0.5 we have the violation of the inequality and 
as ip — > we recover the pure states limit of 0.75. 

In summary, while being admittedly unrealistic, our model offers a few insights into how efficiently the RDDI can 
be used to entangle atoms or implement quantum logic gates (l^Jl^] ■ We have shown that considerable fidelities (up 
to 0.8) of creation of one of the maximally entangled Dicke states and Bell inequality violations can be realized if the 
atoms are placed within distances of the order of a tenth of a wavelength of the working transition. Such distances 
can be achieved, for example, in the ground vibrational states of two atoms in optical lattices pp|Jl^ , [l7|| . 

However, all of the practical applications (say, quantum teleportation [18]]) require stable entanglement. Even the 
Bell inequalities violation considered above can be verified only if the produced entanglement lives sufficiently long so 
that the atoms can be spatially separated for individual addressing and photodetection (imperfectness of the detectors 
constitutes another problem that has not yet been addressed). Of course, this cannot be achieved in our model, since 
the dipole interaction and the decay have the same physical nature and we cannot avoid the latter while making use 
of the former. 

Relatively stable coherences (with lifetimes of the order of seconds) can be generated if we use the Zeeman sublevels 
of the atoms as the working levels (qubits). The RDDI is negligible for them, but using three- level atoms in A- 
configuration instead of two-level atoms we can overcome that difficulty. We can use Raman pulses transferring the 
population between the Zeeman sublevels via a higher lying "transit" Dicke level (STIRAP techniques might be 
an alternative) , so that the transitions between each of the Zeeman sublevels states and the excited "transit" state 
benefit from substantial RDDI. Then by choosing one-photon detunings to be in resonance with only one of the higher 
level Dicke states we can generate the Zeeman sublevels entanglement. 

In this paper we present the first tentative quantitative model of the process of entanglement of atoms with the 
help of the RDDI. Much remains to be done before we can compare the results of the theory with possible practi- 
cal implementations, but the conclusions presented here are promising and therefore encourage further theoretical 
developments. 
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FIG. 1. Geometry of the model with shown directions of laser beams for the symmetric and antisymmetric excitation. 




FIG. 3. a) The fidelity of creation of the maximally entangled states (by applying the optimal pulse as described in text) 
versus interatomic distance ip. b) The l.h.s. of the Bell inequality (5) versus ip. The pure states limit is 0.75, uncorrelated 
states lead to 1.5, the classical boundary is 1.0. Solid and dashed lines correspond to the antisymmetric and symmetric cases, 
respectfully. 
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